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The masses of neutron stars in neutron star binaries are observed to fall in a narrow mass range 
around ~ 1.33 M©. We explore the advantage of focusing on this region of the parameter space in 
gravitational wave searches. We find that an all-sky (externally triggered) search with optimally 
reduced template bank is expected to detect 14% (61%) more binary mergers than without the 
reduction. A reduced template bank can also represent significant improvement in technical cost. 

We also develop a more detailed search method using binary mass distribution, and find similar 
sensitivity increase to that due to the reduced template bank. 


Binary neutron star (NS) mergers represent one of the 
most promising source type for gravitational wave (GW) 
detection [1, 2]. With the recent onset of observations 
with the advanced LIGO detectors [3] and with advanced 
Virgo and KAGRA in the near future [4, 5], the first 
detections are expected within the next few years [1]. 

Gurrent searches for compact binary mergers aim to 
cover virtually the full plausible binary parameter space 
[2] . The strategy is also motivated by the fact that search 
sensitivity is considered to be only weakly dependent on 
the extent of the covered parameter space [6] . 

The primary search method for GWs from compact 
binary mergers is the use of matched filters [2, 6]. This 
method correlates the known signal waveform, called 
template, with the data to identify a GW signal. For 
given binary parameters, the signal waveform can be cal¬ 
culated to high precision. The binary parameter space 
is then covered by using a large number of templates, 
called the template bank, such that the search sensitiv¬ 
ity is sufficiently close to optimal within the considered 
parameter space (e.g., [7]). 

The properties of NS binaries are increasingly con¬ 
strained due to the growing number of observed binaries. 
These observations suggest that the mass of NS within 
NS binaries is within a surprisingly small range around 
~ 1.33 M0, much smaller than the allowed NS mass range 
from ~ 1 Mq to < 3Mq, or even the mass range of NS 
in NS-white dwarf binaries [8]. 

In this paper we investigate the effect of the physically 
motivated reduction of the template bank on search sen¬ 
sitivity. We consider (i) blind, so-called all-sky searches 
(e.g., [9]) in which only GW data is utilized, as well as 
(ii) so-called externally triggered searches (e.g., [10]), in 
which the electromagnetic or other detection of the bi¬ 
nary merger aids the GW search. After discussing the 
dependence of the search sensitivity on the size of the 
template bank, we optimize sensitivity as a function of 
the confidence region of the NS masses in the binary for 
the two search strategies. Finally, we develop a more de¬ 
tailed search that incorporates a ranking statistic for the 
templates based on the expected mass distribution of NS 



FIG. 1. Distribution of GW binary templates in the param¬ 
eter space of the masses of the two NSs in the binary. Gir- 
cles: 92% (smaller) and 99% (bigger) mass confidence regions 
based on the empirical binary NS mass distribution of [8]. 
Templates within the bigger circle (orange+blue) correspond 
to the optimally reduced all-sky search. Templates within 
the smaller circle (blue) correspond to the optimally reduced 
externally triggered search. 


binaries. We calculate the advantage of such search over 
a baseline search that uniformly weights templates. 


Sensitivity dependence on template bank. — Let h{f) 
be the known gravitational waveform and S{{f) be the 
power spectral density of strain noise in the detector. 
The signal-to-noise ratio (SNR) of a matched-filter-based 
search is [6] 
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Assuming the simple case of a single GW detector with 
stationary Gaussian noise, a detection can be claimed 
when the matched filter SNR exceeds a threshold pth- 
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This threshold depends on, among others, the trial factor 
A/trial associated with the search. For Gaussian noise, the 
threshold satisfies [6] 

FAP 

erfc(pth/y2) «-TT-, (2) 

•^vtrial 

where FAP is the false alarm probability of the search. 
LIGO typically considers pth = 8 in a single detector as 
a detection threshold [1]. We adopt this value in the 
following. 

The trial factor A/triai depends on the number 
Atempiates of templates used, along with the number Aft 
of independent starting times that are considered [6] , and 
for analyses with multiple detectors, the number A/n of 
independent source directions the search considers. Us¬ 
ing the fact that At oc fobs, where fobs is the observation 
duration, and A/n oc fl, where is the allowed sky region, 
the trial factor satisfies 


—0.03. We adopt this empirical NS mass distribution in 
the following. We conservatively assume that the two NS 
masses in the binary are independent. 

We utilize the NS mass probability distribution by us¬ 
ing only those templates in a matched-filter-based search 
that are more likely to be observed from astrophysical 
sources. This way, we can significantly reduce the size 
of the template bank with modest reduction of the frac¬ 
tion /ns of NS binaries whose masses are covered by the 
template bank. Let Yns be the region in the 2D mass 
parameter space that is included in the analysis. For any 
/ns, we can define Uns such that (i) a template that is 
S F/is is at least as likely to correspond to a detected NS 
binary as any template that is ^ Fns, and (ii) the fraction 
of NS binaries that fall within Fns is /ns- 

Sensitivity improvement. — The most meaningful 
quantity to compare searches with is the detection rate 
TZ, which can be written as 


A/trial CC A/templates ^obs (3) 

We see that changing the size of the template bank 
changes A/triai, which in turn will change the search sen¬ 
sitivity through Pth- 

Besides the effect of the template bank, we also see 
that fobs and D similarly affect search sensitivity, which 
is relevant for externally triggered searches. 

We now derive a formula to quantify the change in 
sensitivity due to varying the search parameters. Let our 
baseline search have -A^te^piates templates, observa¬ 

tion time, include sky area, and have po detection 
threshold. We want to determine the new threshold pi 
if we change the parameters to -/V^empiates’ ^ibi 
respectively. We can use Eq. 2, which yields, to a good 
approximation, pth ~ y^2 ln(Atriai/FAP) [6]. With this, 
we obtain 
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utilizing neutron star mass distribution. — The 
masses of NSs within observed NS binaries fall in a sur¬ 
prisingly small range. While the maximum allowed NS 
mass is above 2Mq [12], NS masses in NS binaries are 
closely clustered around 1.33 Mq. Kiziltan et al. [8] use 
the observed masses and their uncertainties in a statisti¬ 
cal model to find an empirical NS mass (m) distribution 


= (5) 

where (/(a;) and $(a;) are the standard normal density 
and cumulative density functions, respectively. For NS 
binaries, Kiziltan et al. hnd p = 1.33, cr = 0.11 and a = 


TZ — fnsfh /gw- (6) 

Here, ft is the beaming factor of the emission correspond¬ 
ing to the external trigger, which reduces the number of 
observed sources [13]; for all-sky searches, we can take 
/{, = !. For externally triggered searches, the direction 
of the external trigger can be correlated with the weakly 
direction-dependent GW emission. For triggers such as 
GRBs, which are aligned with the orbital axis of the bi¬ 
nary merger, for /{, ^ 1 we can approximate /gw « 1.5 
[13], which is the ratio of the GW strain amplitude in the 
direction of the orbital axis compared to the directionally 
averaged strain amplitude. 

To estimate the sensitivity improvement in different 
search scenarios, we adopt the template bank used by 
initial LIGO-Virgo in [9], which is representative of what 
will be used for advanced LlGO-Virgo searches [2] (see 
Fig. 1 for the distribution of templates for the mass 
regime < 2 M©). The templates are distributed such that 
the loss in SNR due to the discreteness of the template 
bank is less than 3%, but otherwise cover the parameter 
space for masses > 1 Mq . To demonstrate the role of 
the size of the template bank for binary NS searches, we 
focus on the part of the bank for which both masses in the 
binary are < 3 Mq . This is the parameter range planned 
to be used for binary NS searches with advanced LIGO 
[2] . The part of the template bank below 3 Mq includes 
~ 45,000 templates. We denote the detection rate of this 
baseline search with 7?,aiisky- 

We first consider all-sky observations. We assume 
fobs = 1 yr observation time, and O = 47r. The baseline 
search, i.e. the one without taking advantage of the NS 
mass distribution, will use the full template bank below 
3M0. We take a detection threshold of Paiisky = 8 [1]. 
Compared to this baseline, an all-sky search using a con¬ 
strained binary NS parameter space will yield detection 
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FIG. 2. Improvement of detection rate using the constrained 
NS mass parameter space over the baseline search with no 
constraints.Results are shown for all-sky (solid) and externally 
triggered (dashed) searches, as functions of the fraction /ns of 
the NS mass parameter space included in the template bank. 


threshold Paiisky.NS < Paiisky due to the reduced template 
bank. The detection rate for the constrained parameter 
space, normalized by the detection rate for the baseline 
all-sky search, is shown in Fig. 2 as a function of the NS 
confidence level /ns- We select the /ns value that max¬ 
imizes the detection rate 7?.aUsky,NS- The corresponding 
value is /ns.aiisky = 0.99, with Paiisky.NS = 7.6 and rate 
TT-aiisky.Ns/TT-aiisky = 1.14. The Constrained parameter 
space and the corresponding template distribution are 
shown in Fig. 1. 

We next look at externally triggered searches.External 
triggers can significantly boost search sensitivity by de¬ 
creasing the time window and sky area in which GWs 
need to be searched for. Similarly to Chen and Holz [11], 
we consider fobs = 10 s, corresponding to the approxi¬ 
mate GW search time window for one sufficiently nearby 
GRB detected, and ft = 100 deg^, an approximate sky 
area corresponding to a known signal direction, taking 
into account the directional precision of GW searches. 
For the baseline externally triggered search, we do not 
take advantage of the NS mass distribution, and use the 
full template bank below SM©. Using Eq. 4, we find 
that the detection threshold of this baseline search is 
Pextrig = 4.7. We now obtain the detection rate improve¬ 
ment for the case with constrained NS parameter space. 
Similarly to the all-sky case, we select the /ns value that 
maximizes the detection rate 77.extrig,NS- The correspond¬ 
ing value is /ns.extrig = 0.92, with Pextrig.NS = 3.9 and rate 
77-oxtrig,Ns/77-oxtrig = 1.61, a significant improvement over 
the baseline externally triggered case. The constrained 
parameter space and the corresponding template distri¬ 
bution are shown in Fig. 1. 

We now explore the difference between detection rates 
for the all-sky and externally triggered cases. This 
difference depends on the beaming factor ft = [1 — 


cos(0)]“^, where 9 is the opening angle. In general, ex¬ 
ternally triggered searches benefit from increased sensi¬ 
tivity, while all-sky searches are advantageous as only a 
fraction of events will have external triggers due to beam¬ 
ing. Using Eq. 6, we find that, for the baseline case, 
TT-extrig/TT-aiisky = 0.25 and 2.24 for 9 = 10° and 30°, re¬ 
spectively. The two baseline searches are expected to de¬ 
tect the same number of binaries for 9 = 20°. For the con¬ 
strained searches, TT-extrig.Ns/TT-aiisky.NS = 0.36 and 3.18 
for 9 = 10° and 30°, respectively, significantly higher 
than for the baseline case. The constrained externally 
triggered search is expected to do as well as the con¬ 
strained all-sky search for 9 = 17°. This means that 
for the typically expected short-GRB opening angles of 
< 10° [14], the majority of the detected events will come 
from all-sky searches. 

Search method incorporating NS mass distribution. — 

To gain a more detailed picture of the advantages of 
utilizing the binary NS mass distribution in GW searches, 
we calculate the increased sensitivity of a likelihood-ratio 
test that incorporates mass information (for likelihood- 
based methods and their benefits, see [15-20]). We define 
our test statistic as 

^ _ P(toi, m 2 |signal)P(p|TOi,m 2 ,signal) 

P(toi, m2|noise)P(/5|TOi, m 2 , noise) 

where mi and m 2 are the masses of the two NS in the bi¬ 
nary, and P(-1 signal/noise) is the conditional probability 
for the signal/noise hypothesis. Here, P(mi,m 2 |-) is the 
probability corresponding to the template with masses 
mi and m 2 , i.e. it is not a probability density. 

For simplicity, we assume P(mi, m 2 [noise) = 
const, and that the probability distribution of p is 
mass independent: P(p|mi,m 2 ,noise) = P(/o|noise) 
and P(p|mi,m 2 ,signal) = P(p|signal). Note that 
P(mi,m 2 |-) considered here is the probability of a spe¬ 
cific template defined by mi and m 2 , i.e. not probability 
density. 

Adopting a uniform binary NS distribution in the local 
universe yields P(/9|signal) oc p~'^. Gonsidering a Gaus¬ 
sian noise model, we get P(p|noise) oc exp(—p^/2). With 
these assumptions, we arrive at 

£ = P(mi, m 2 |signal)(8) 

where we omitted a constant factor that is irrelevant for 
the statistic. 

We consider three different models for 
P(mi, m 2 [signal). The baseline model that as¬ 
signs no weight to the specific {mi, m 2 } val¬ 
ues: P(mi, m 2 [signal) [baseline = coTist. A model 

with a cutoff in the allowed mass parameter 
space, in line with our discussion above; for 

this model, P(mi, m 2 [signal) [cutoff = const, for 
+ — < 3cr and = 0 otherwise. In the 





























4 


third, weighted model, we fully make use of the expected 
mass distribution by assigning the appropriate weight 
fromEq. 5: P(mi, m 2 |signal)|weighted =-Pns(wi)Pns(TO 2 ) 
(see e.g., [18] for context and motivation of weighted 
signal priors). 

We are interested in calculating the sensitivity of the 
different models as functions of FAP. For this, we define 
FAP as the probability of obtaining £ greater than a 
threshold value £th during the observation period for any 
of the templates: FAP = Pobs('^^ > £th[noise). Instead 
of £th, we can also equivalently specify the threshold as 
a SNR threshold pth using Fq. 8. 

/ OO _ 

P(/9|noise)dp = Eli-i-tai...)] (9) 

■i rth,Z 2 

where pth,i is the threshold value of the template such 
that £(mi_i, m 2 y, pth,i) = £th, the sum is over all tem¬ 
plates in the bank, and <i> is the cumulative distribution 
function of the standard normal distribution. 

To calculate the search sensitivity as a function of 
FAP, we calculate the expected rate of detected NS bi¬ 
naries. We consider a binary detected if its SNR is 
> /Oth(wi, TO 2 ), where we note that the threshold depends 
on the binary masses. For uniform source distribution, 
the expected rate of detected binaries is oc p~^. For sim¬ 
plicity, we neglect the effect of decreased sensitivity due 
to binary masses being off of templates, and consider that 
the template with the closest masses to those of the bi¬ 
nary will be used. We define distance with respect to 
mass with (Amf -I- Am|)^/^. Using this definition, for 
template i we determine the area Ai within the mass pa¬ 
rameter space that is closest to template i. Assuming 
that, for a given binary merger, its SNR and the closest 
template can be accurate recovered, the total detected 
rate can be expressed as 

7^oc m 2 ,i[signal) A, (10) 

i 

Using Eqs. 9 and 10, we can calculate the expected detec¬ 
tion rate as a function of FAP for each of the considered 
search models. The results are shown in Fig. 3. 

To compare the sensitivity of different search mod¬ 
els, we select FAPq that corresponds to pth = 8 for the 
baseline all-sky search, the typical value used by LIGO 
searches [1], and the one considered above for the effect 
of template bank reduction. Note that pth is identical 
for all templates for the baseline search. We calculate 
the corresponding pth for externally triggered searches 
similarly to the analysis of the effect of reduced template 
banks above. The results in Fig. 3 are shown such that 
FAP values are normalized by FAPq, and the expected 
detection rates are normalized by the detection rates TZq 
corresponding to those of the baseline searches. This al¬ 
lows for the quantification of the advantage of utilizing 
the expected NS mass distribution in the test statistic. 



FIG. 3. Expected detection rate as a function of the false 
alarm rate for different all-sky and externally triggered search 
models using the likelihood-ratio test statistic, normalized to 
the baseline method. 

The ’cutoff’ model comparison represents the same 
improvement as the reduced-template-bank comparison 
above. We find that the expected detection rate for the 
cutoff model improves by ^ 12% and 35% for the all-sky 
and externally triggered cases compared to the baseline, 
respectively. For the ’weighted’ model, we see improve¬ 
ments of ^ 17% and 68% for the two cases, respectively. 

Conclusion. — We explored the increase in the de¬ 
tection rate of GWs from binary NS mergers due to re¬ 
ducing the search template bank based on the expected 
mass distribution of NS in binaries. We found an increase 
of 14% for all-sky searches, and 61% for an externally 
triggered search. The higher increase for the externally 
triggered case was expected since in that case the tem¬ 
plate bank is a more significant contributor to the search 
trial factor. For both all-sky and externally triggered 
cases, the optimal detection strategy includes templates 
that cover over 90% of the expected binary NS mergers, 
therefore increasing the detection rate can be achieved 
by covering a significant fraction of the binary NS mass 
parameter space, i.e. hardly losing any astrophysical sig¬ 
nal. 

We constructed realistic search method that incorpo¬ 
rates the expected binary NS mass distribution in its test 
statistic, and calculated its sensitivity increase over using 
no information about the mass distribution. We find that 
this realistic search achieves similar sensitivity increases 
for both all-sky and externally triggered searches than 
what we found by simply reducing the template bank. 

Beyond the obvious benefit of increasing the detection 
rate, constraining the NS template bank also significantly 
reduces the computational cost of the search, with the 
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required template bank being only 6% and 3% of the total 
template bank for the all-sky and externally triggered 
cases, respectively. 

Several future extensions of the present work will be 
interesting. Further observations of NS binary masses 
can help, e.g. through establishing correlations between 
the two NS masses. While current masses are for nearby 
binaries, gaining farther sample can improve detection 
prospects [17]. Additionally, similar improvements can 
improve black hole-NS and black hole-black hole mergers 
as well. 

The authors thank Kipp Cannon, Tom Dent, John 
Veitch and Zsuzsa Marka for their useful suggestions. 
This paper has been approved for publication by the 
LIGO Scientific Collaboration, with LIGO Document 
Number LIGO-P1500134. IB and SM are thankful for 
the generous support of Columbia University in the City 
of New York and the National Science Foundation under 
cooperative agreement PHY-1404462. 


* ibartos@phys.columbia.edu 

[1] LIGO Scientific Collaboration et ai, ArXiv e-prints 
(2013), 1304.0670. 

[2] LIGO Scientific Collaboration and Virgo Collaboration, 
The LSC-Virgo White Paper on Gravitational Wave 
Searches and Astrophysics, 2014, https: //dec. ligo. 
org/LIG0-T1400054/public. 

[3] The LIGO Scientific Collaboration et al., Classical and 


Quantum Gravity 32, 074001 (2015), 1411.4547. 

[4] F. Acernese et al, Classical and Quantum Gravity 32, 
024001 (2015), 1408.3978. 

[5] The KAGRA Collaboration, Y. Aso et al, Phys. Rev. D 
88, 043007 (2013). 

[6] E. E. Flanagan and S. A. Hughes, Phys. Rev. D57, 4535 
(1998), gr-qc/9701039. 

[7] D. A. Brown, 1. Harry, A. Lundgren, and A. H. Nitz, 
Phys. Rev. D86, 084017 (2012), 1207.6406. 

[8] B. Kiziltan, A. Kottas, M. De Yoreo, and S. E. Thorsett, 
Astrophys. J.778, 66 (2013), 1011.4291. 

[9] J. Abadie et al., Phys. Rev. D85, 082002 (2012), 
1111.7314. 

[10] B. Abbott et al, Astrophys. J.681, 1419 (2008), 
0711.1163. 

[11] H.-Y. Chen and D. E. Holz, Physical Review Letters 
111, 181101 (2013), 1206.0703. 

[12] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. 
Roberts, and J. W. T. Hessels, Nature (London)467, 
1081 (2010), 1010.5788. 

[13] 1. Bartos, P. Brady, and S. Marka, Classical and Quan¬ 
tum Gravity 30, 123001 (2013), 1212.2289. 

[14] W. Pong et al, Astrophys. J.756, 189 (2012), 1204.5475. 

[15] R. Biswas et al., Phys. Rev. D85, 122008 (2012), 
1201.2959. 

[16] R. Biswas et al., Phys. Rev. D85, 122009 (2012), 
1201.2964. 

[17] T. Dent and J. Veitch, Phys. Rev. D89, 062002 (2014), 
1311.7174. 

[18] K. Cannon, C. Hanna, and J. Peoples, ArXiv e-prints 
(2015), 1504.04632. 

[19] J. Abadie et al., Phys. Rev. D82, 102001 (2010), 
1005.4655. 

[20] J. Abadie et al., Phys. Rev. D85, 089903 (2012). 



